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We study nonequilibrium phase transitions in a mass-aggregation model which allows for diffu- 
sion, aggregation on contact, dissociation, adsorption and desorption of unit masses. We analyse 
two limits explicitly. In the first case mass is locally conserved whereas in the second case local 
conservation is violated. In both cases the system undergoes a dynamical phase transition in all 
dimensions. In the first case, the steady state mass distribution decays exponentially for large mass 
in one phase, and develops an infinite aggregate in addition to a power-law mass decay in the other 
phase. In the second case, the transition is similar except that the infinite aggregate is missing. 
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The steady state of a system in thermal equilibrium 
is described by the Gibbs distribution. Phase transi- 
tions which occur in such equilibrium systems as one 
changes the external fields such as temperature or mag- 
netic field are by now weU understood. On the other 
hand there is a wide variety of inherently nonequilib- 
rium systems in nature whose steady states are not de- 
scribed by the Gibbs distribution, but are determined 
by the underlying microscopic dynamical processes and 
are often hard to determine. Examples include systems 
exhibiting self-organized criticality several reaction- 
diffusion systems fluctuating interfaces |^ and many 
others. As one changes the rates of the underlying dy- 
namical processes, the steady states of such systems may 
undergo nonequilibrium phase transitions. As compared 
to their equilibrium counterparts, these nonequilibrium 
steady states and the transitions between them are much 
less understood due to the lack of a general framework. 
It is therefore important and necessary to study simple 
models amenable to analysis in order to understand the 
mechansims of such phase transitions. 

Here we study the nonequilibrium phase transitions in 
an important class of systems which involve microscopic 
processes of diffusion and aggregation, dissociation, ad- 
sorption and desorption of masses. These processes are 
ubiquitous in nature, and arise in a variety of physical 
settings, for example, in the formation of colloidal sus- 
pensions [Q and polymer gels on the one hand, and 
aerosols and clouds on the other. They also enter in an 
important way in surface growth phenomena involving is- 
land formation In this Letter, we introduce a simple 
lattice model incorporating these microscopic processes 
and study the nonequilibrium steady states and the tran- 
sitions between them both analytically within mean field 
theory and numerically in one dimension. 

Our lattice model, which evolves in continuous time, is 
defined as follows. For simplicity we define the model on 
a one-dimensional lattice with periodic boundary condi- 
tions although generalizations to higher dimensions are 
quite straightforward. Beginning with a state in which 
the masses are placed randomly, a site is chosen at ran- 



dom. Then one of the following events can occur: 

1. Adsorption: With rate q, a unit mass is adsorbed 
at site i; thus mi nii + 1. 

2. Desorption: With rate p, a unit mass desorbs from 
site i: thus TOj — )■ to^ — 1 provided to^ > 1. 

3. Chipping (single-particle dissociation): With rate 
w, a bit of the mass at the site "chips" off, i.e. 
provided rui > 1, a single particle leaves site i and 
moves with equal probability to one of the neigh- 
bouring sites i — I or i + 1 ; thus — > — 1 and 
mj±i -I- 1. 

4. Diffusion and Aggregation: With rate 1, the mass 
mi at site i moves either to site z — 1 or to site i + 1. 
If it moves to a site which already has some parti- 
cles, then the total mass just adds up; thus — ^ 
and mi±i mi±i + mi. 

Note that we have assumed that both desorption and 
diffusion rates are independent of the mass. In a more re- 
alistic situation these rates would depend upon the mass. 
However, our aim here is not to study this model in full 
generality, but rather to identify the mechanism of a dy- 
namical phase transition in the simplest possible scenario 
involving these microscopic processes. Indeed we show 
below that even within this simplest scenario, novel dy- 
namical phase transitions occur which are nontrivial yet 
amenable to analysis. 

Though the model can be studied in the full parameter 
space of all four basic processes, for simplicity we restrict 
ourselves here to two limiting cases: (i) p = 0, g = 0, 
i.e. only chipping, diffusion and aggregation moves are 
allowed. In this limit, mass is locally conserved by the 
moves and we call this model the conserved-mass aggre- 
gation model (CMAM) (ii) w = 0, i.e. all moves except 
for chipping are allowed. In this case, adsorption and 
desorption lead to violation of local mass conservation. 
We call this the In-out model. In this Letter, we analyse 
the CMAM model in some detail and only outline the 
main results for the In-out model. 
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Let us summarize our main results: (i) In the CMAM, 
single particles are allowed to chip off from massive con- 
glomerates. This move corresponds to the physical pro- 
cess of single functional units breaking off from larger 
clusters in the polymerization problem. It leads to a re- 
plenishment of the lower end of the mass spectrum, and 
competes with the tendency of the coalescence process 
to produce more massive aggregates. The result of this 
competition is that two types of steady states are possi- 
ble, and there is a dynamical phase transition between 
the two. In one state, the steady state mass distribu- 
tion P{in) decays exponentially, while the other is more 
striking and interesting: P{m) decays as a power law for 
large m but in addition develops a delta function peak at 
TO = oo. Physically this means that an infinite aggregate 
forms that subsumes a finite fraction of the total mass, 
and coexists with smaller finite clusters whose mass dis- 
tribution has a power law tail. In the language of sol-gel 
transitions, the infinite aggregate is like the gel while the 
smaller clusters form the sol. However, as opposed to the 
models of irreversible gelation where the sol disappears 
in the steady state, in our model the sol coexists with the 
gel even in the steady state. Interestingly, the mechanism 
of formation of the infinite aggregate in the steady state 
resembles Bose-Einstein condensation (BEC), though the 
condensate (the infinite aggregate here) forms in real 
space rather than momentum space as in conventional 
BEC. (ii) In the In-out model too we find a phase tran- 
sition in the steady state as the adsorption (q) and des- 
orption (p) rates are varied. In one phase (low values of 
q) P{m) decays exponentially whereas in the other phase 
(high q) it has a power law tail. This power law phase 
is similar to that of the Takayasu model ||^ of particle 
injection and aggregation. 

We first analyse the CMAM within the mean field 
approximation, ignoring correlations in the occupancy of 
adjacent sites. Then we can directly write down equa- 
tions for P{m, t), the probability that any site has a mass 
TO at time t. 



dP{m,t) 
dt 

dP{0,t) 
Jt 



= -(1 + + s{t)]P{m,t) + wP{m + l,t) 



wsit)P{m-l,t) + P*P; TO>1 



(1) 



-(1 + w)sit)P{0, t) + wP{l, t) + s{t). (2) 



Here s{t) = 1 — P(0, t) is the probability that a site is oc- 
cupied by a mass and P*P = X]m'=i P{^', t)P{m—m', t) 
is a convolution term that describes the coalescence of 
two masses. 

The above equations enumerate all possible ways in 
which the mass at a site might change. The first term in 
Eq. (|l|) is the "loss" term that accounts for the probabil- 
ity that a mass to might move as a whole or chip off to 
either of the neighbouring sites, or a mass from the neigh- 
bouring site might move or chip off to the site in consider- 



ation. The probability of occupation of the neighbouring 
site, s{t) = X]m=i ^('^' multiplies P{m,t) within the 
mean-field approximation where one neglects the spatial 
correlations in the occupation probabilities of neighbour- 
ing sites. The remaining three terms in Eq. (^ arc the 
"gain" terms enumerating the number of ways that a site 
with mass m' ^ to can gain the deficit mass to — to'. The 
second equation Eq. (|2|) is a similar enumeration of the 
possibilities for loss and gain of empty sites. Evidently, 
the mean field equations conserve the total mass. 

To solve the equations, we compute the generating 
function, Qiz,t) = Em=i *)^'" ^om Eq. (0) and 
set dQ/dt = in the steady state. We also need to use 
Eq. (||) to write P(l,i) in terms of s{t). This gives us a 
quadratic equation for Q in the steady state. Choosing 
the root that corresponds to Q(z = 0) = 0, we find 



w 



2s 



ws- 



(1 



2z 



ws w wsz 
2l 2~ 

■y/iz - zi){z - Z2). 



(3) 



where zi^2 = (w -I- 2 =F 2y/w + l)/ws. The value of the 
occupation probability s is fixed by mass conservation 
which implies that ^ mP{m) — M/L = p. Putting 
dzQ{z = 1) = p, the resulting relation between p and 
s is 



2p = w{l - s) - WSy/ (zi - 1)(Z2 - 1) . 



(4) 



The steady state probability distribution P{m) is the 
coefficient of z™ in Q{z) and can be obtained from Q{z) 
in Eq. (^ by evaluating the integral 



P(m) 



1 
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dz 



(5) 



over the contour Co encircling the origin. The singular- 
ities of the integrand govern the asymptotic behaviour 
of P{m) for large to. Clearly the integrand has branch 
cuts at z = zi_2. For fixed w, if one increases the density 
p, the occupation probability s also increases as evident 
from Eq. (||). As a result, both the roots zi^2 start de- 
creasing. As long as the lower root Zi is greater than 1, 
Eq. ^ is well defined and the analysis of the contour 
integration around the branch cut z = zi , yields for large 

TO, 



P(to) 



-m/m* ^^3/2 



(6) 



where the characteristic mass, to* = l/log(zi) and di- 
verges as ^ (sc — s)~^ as s approaches Sc = {w + 2 — 
2\/w + 1)/^. Sc is the critical value of s at which zi = 1. 
This exponentially decaying mass distribution is the sig- 
nature of "disordered" phase which occurs for s < Sc or 
equivalently from Eq. (0) for p < pc = y/w + 1 — 1. 
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When p = pc, we have zi = 1, and analysis of the 
contour around z — zi — 1 yields a power law decay of 
P{m), 



P(m) 



-5/2 



(7) 



As p is increased further beyond pc, s cannot increase 
any more because if it does so, the root zi would be less 
than 1 (while the other root Z2 is still bigger than 1) and 
Eq. would be undefined. The only possibility is that s 
sticks to its critical value Sc or equivalently the lower root 
zi sticks to 1. Physically this implies that adding more 
particles does not change the occupation probability of 
sites. This can happen only if all the additional particles 
(as p is increased) aggregate on a vanishing fraction of 
sites, thus not contributing to the occupation of the oth- 
ers. Hence in this "infinite-aggregate" phase P{m) has 
an infinite-mass aggregate, in addition to the power law 
decay vnT^I"^ . Concomitantly Eq. (^) ceases to hold, and 
the relation now becomes 



(8) 



where poo is the fraction of the mass in the infinite ag- 
gregate. The mechanism of formation of the aggregate is 
reminiscent of Bose Einstein condensation. In that case, 
for temperatures in which a macroscopic condensate ex- 
ists, particles added to the system do not contribute to 
the occupation of the excited states; they only add to the 
condensate, as they do to the infinite aggregate here. 

Thus the mean field phase diagram (see inset of Fig. 1) 
of the system consists of two phases, "Exponential" and 
"Aggregate", which are separated by the phase bound- 
ary, Pc = \fw + \ — 1 . While this phase diagram remains 
qualitatively the same even in 1-d, the exponents char- 
acterizing the power laws are different from their mean 
field values (see Fig. 1). 
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1. log-log plot of P(m) vs. m for the CMAM, for 
> Pc on a lattice with L — 1000. Inset: Phase dia- 
gram. The solid line and the points indicate the phase bound- 
ary within mean field theory and 1-d simulation respectively. 
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We have studied this model using Monte Carlo simula- 
tions on a one-dimensional lattice. Although we present 
results here for a relatively small size lattice, L = 1024, 
we have checked our results for larger sizes as well. 
We confirmed that all the qualitative predictions of the 
mean-field theory remain true, by calculating P(m) nu- 
merically in the steady state. Figure 1 displays two nu- 
merically obtained plots of P{m). The existence of both 
the "Exponential" (denoted by x) for p < pc and the 
"Aggregate" phase (denoted by for p > pc is con- 
firmed. In particular, the second curve shows clear evi- 
dence of a power-law behaviour of the distribution, which 
is cut off by finite-size effects, and for an 'infinite' ag- 
gregate beyond. We confirmed that the mass Magg in 
this aggregate grows linearly with the size, and that the 
spread 6Magg grows sublinearly, implying that the ra- 
tio S Magg /Magg approachcs zero in the thermodynamic 
limit. The exponent t^mam which characterizes the 
finite-mass fragment power law decay for p > pc is nu- 
merically found to be 2.33 ± .02 and remains the same at 
the critical point p = pc- 

We note that in conserved-aggregation models studied 
earlier within mean field theory [^,0, the steady state 
mass distribution also changed from an exponential dis- 
tribution to a power law as the density was increased 
beyond a critical value. However, the existence of the 
striking infinite aggregate in the steady state for p > Pc 
was not identified earlier. 

We next study the steady state phase diagram of the 
In-out model in the q-p plane. In this model, mass is 
not locally conserved. The mass per site M{t) evidently 
obeys the exact equation 



dM 



(9) 



where s{t) is the fraction of sites occupied by a mass 
mi > 1. In the steady state, let the mean value of s{t) 
be s. If q/p is low, s adjusts to make q — ps vanish, 
and the mean mass reaches a time- independent value M. 
This defines the finite-mass phase. As we will see below, 
as q/p increases beyond a critical value, s never catches 
up with q/p and reaches a steady state value which is 
less than q/p; in this phase, M increases linearly in time 



while P{m, t) 



f{mt ^) which in the long time 



limit converges to a time-independent form, decaying as 
a power law with exponent tt, even though the moments 
of this distribution diverge as time increases to infinity. 
We call this the growing-mass or the Takayasu phase. In 
fact, for p — the In-out model reduces exactly to that of 
the Takayasu model (TM) of injection and aggregation 
of masses [|j which has found widespread applications 
ranging from river models to granular systems p2[ . 
Indeed what we find here is that the growing mass phase 
of the TM at p = persists up to a nonzero critical value 
Pciq) for a given q, while for p > Pc{q) the mass stops 
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growing and P{m) decays exponentially for large m in 
the steady state. 

The mean field analysis of the In-out model is similar 
to that of the CMAM model though a little bit trick- 
ier. We defer the details for a future publication and 
only outline the results here. We find that the critical line 
Pc{q) = 1 + 2^/? separates two phases in the q-p plane. 
For p > pc, P{m) ^ mr'^l'^ exp(— m/m*) for large m. For 
p = Pc, P{m) ^ m~^/^ and for p < pc, P{m) ~ mT^/"^ 
for large m. For a fixed q, the steady state occupation 
density s{p,q) develops an interesting cusp as p crosses 
Pc{q)- For example, at g = 1, where pc = 3, s{p) = 1/p 
for p > 3 as follows simply from Eq. (9); but for p < 3, 
the determination of s{p) is nontrivial and is given 
by the positive root of the cubic equation, 16ps^ + (8p^ + 
4p-25)s2 + (p^-llp2-43p-25)s-p3+2p2 + i7p+25 = 0. 

The qualitative predictions of mean field theory - the 
existence of a power-law (Takayasu) phase (P(m) ~ 
mT'^'^) and a phase with exponential mass distribu- 
tion, with a different critical behaviour at the transition 
{P(m) ^ 771"'^'=) - are found to hold in 1-d as well. The 
Takayasu exponent tt is known exactly to be 4/3 in 1-d 
and 3/2 within mean field theory Figure 2 shows the 
results of numerical simulations in 1-d for the phase dia- 
gram and the decay of the mass distribution in the two 
phases and at the transition point. The values obtained, 
Tt = 4/3 and Tc ~ 1.833, are quite different from their 
mean-field values, tt — 3/2 and Tc = 5/2, reflecting the 
effects of correlations between masses at different sites. 
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FIG. 2. log-log plot of P{m) vs. m for the In-out model, 
for p <,—,> Pc- Inset: Phase diagram. The solid line and the 
points indicate the phase boundary within mean field theory 
and 1-d simulation respectively. 

We may reinterpret the configuration of masses in the 
In-out model as an interface profile on regarding rrii as 
a local height variable. While the model may have some 
unphysical features in the context of an interface due to 



the columns of masses moving as a whole, the analogy 
helps however to understand physically the nature of the 
transition in the In-out model. In the interface language 
this corresponds to a wetting transition; the key factor 
responsible for the occurrence of the smooth phase is a 
substrate, implicit in the constraint rTi^ > in the In- 
out model. The wet phase is identified with a growing 
mass phase, which has a rough profile, with exact rough- 
ness exponent xt — 5/2 Q in 1-d. Since xt > 1, the 
interface in the wet phase is not self-affine. Recently a 
nonequilibrium wetting transition was also observed in 
an interface model where the interface in the wet 
phase is self-affine due to surface tension effects which 
are absent in our model. Interestingly, however, in our 
model the substrate is able to induce a self-affine interface 
at the critical point with roughness exponent x — 1/3 
within mean field theory and x ~ 0.7 in 1-d |p^ , despite 
the anomalously large roughness of the wet phase. 

We thank Deepak Dhar for useful discussions and S. 
Cueille and S. Rcdner for pointing out references |^ and 
to us. 
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